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We study the phases and phase transition hnes of the finite temperature G2 Higgs model. Our work is based on 
an efficient local hybrid Monte-Carlo algorithm which allows for accurate measurements of expectation values, 
histograms and susceptibilities. On smaller lattices we calculate the phase diagram in terms of the inverse gauge 
coupling /3 and the hopping parameter n. For k — ^ the model reduces to G2 gluodynamics and for k — > 00 to 
517(3) gluodynamics. In both limits the system shows a first order confinement-deconfinement transition. We 
show that the first order transitions at asymptotic values of the hopping parameter are almost joined by a line 
of first order transitions. A careful analysis reveals that there exists a small gap in the line where the first order 
transitions turn into continuous transitions or a cross-over region. For /? — > 00 the gauge degrees of freedom 
are frozen and one finds a nonlinear 0(7) sigma model which exhibits a second order transition from a massive 
0(7)-symmetric to a massless 0(6)-symmetric phase. The corresponding second order line for large /3 remains 
second order for intermediate /3 until it comes close to the gap between the two first order lines. Besides this 
second order line and the first order confinement-deconfinement transitions we find a line of monopole-driven 
bulk transitions which do not interfer with the confinement-deconfinment transitions. 

PACS numbers: 11.15.-q, ll.15.Ha, 12.38.Aw 



Quarks and gluons are confined in mesons and baryons and are not seen as asymptotic states of strong interaction. Understanding 
the dynamics of this confinement mechanism is one of the challenging problems in strongly coupled gauge theories. Confinement 
is lost under extreme conditions: when temperature reaches the QCD energy scale or the density rises to the point where the 
average inter-quark separation is less than 1 fm, then hadrons are melted into their constituent quarks. 
For gauge groups with a non-trivial center is the Polyakov loop 



an order parameter for the transition from the confined to the unconfined phase in gluodynamics (pure gauge theories). Its 
thermal expectation value is related to the difference in free energy F due to the presence of an infinitely heavy test quark in the 
gluonic bath as 



such that (P) 7^ in the unconfined high-temperature phase and (P) = in the confined low-temperature phase. Below 
the critical temperature is P{x) uniformly distributed over the group manifold and above the critical temperature it is in the 
neighborhood of a center-element. Near the transition point its dynamics is successfully described by effective three dimensional 
scalar field models for the characters of V (x) [1-3] . If one further projects the Polyakov loops onto the center of the gauge group, 
then one arrives at generalized Potts models describing the effective Polyakov-loop dynamics [4]. 

With matter in the fundamental representation the center symmetry is explicitly broken and for all temperatures has P a non- 
zero expectation value and points in the direction of a particular center element. Thus in the strict sense the Polyakov loop 
ceases to be an order parameter for the center symmetry. On a microscopic scale this is attributed to the breaking of the string 
connecting a static 'quark anti-quark pair' when one tries to separate the static charges [5]. It breaks via the spontaneous creation 
of dynamical quark anti-quark pairs which in turn screen the individual static charges. 

To clarify the relevance of the center symmetry for confinement it suggests itself to study gauge theories for which the gauge 
group has a trivial center Then the Polyakov loop ceases to be an order parameter even in the absence of dynamical matter 
since the strings connecting external charges can break via the spontaneous creation of dynamical 'gluons'. The smallest simple 
and simply connected Lie group with a trivial center is the 14 dimensional exceptional Lie group G2. This is one reason why 
G2 gauge theory with and without Higgsfields has been investigated in series of papers [6-11]. Although there is no symmetry 
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reason for a deconfinement phase transition in G2 gluodynamics it has been conjectured that a first order deconfinement transition 
without order parameter exists. In this context confinement refers to confinement at intermediate scales, where a Casimir scahng 
of string tensions has been detected in [12]. Although the threshold energy for string breaking in G2 gauge theory is rather high, 
string breaking has been seen in 3 dimensional G2 gluodynamics in [13]. 

The gauge group SU (3) of strong interaction is a subgroup of G2 and this observation has interesting consequences, as pointed 
out in [8]. With a Higgs field in the fundamental 7 dimensional representation one can break the G2 gauge symmetry to the 
SU{3) symmetry via the Higgs mechanism. When the Higgs field in the action 

S[A,(t,]= J d*x(^^tTF^,F^^'' + ^{D^c^,D^c^) + V{c^)^ (3) 

picks up a vacuum expectation value v, then 6 gauge bosons acquire a mass proportional to v while the 8 gluons belonging to 
SU{3) remain massless. The massive gauge bosons are removed from the spectrum for v 00. In this limit G2 Higgs model 
reduces to SU{3) Yang-Mills theory. Even more interesting, for intermediate and large values of v the G2 Yang-Mills-Higgs 
(YMH) theory mimics SU{3) gauge theory with dynamical 'scalar quarks'. The masses of these 'quarks' and the length scale at 
which string breaking occurs increase with increasing v. The Polyakov loop serves as approximate order parameter separating 
the confined from the unconfined phases with a rapid change at the transition or crossover This rapid change is depicted in 
Fig. 1 which shows the expectation value of P for G2 gluodynamics as function of the inverse gauge coupling j3 = 1/ g^. 




FIG. 1: Phase transition on a 16^ x 6 lattice in terms of the Polyakov loop in the fundamental representation of G2. The rapid change of the 
Polyakov loop with /3 — (left panel) and histograms of the Polyakov loop for /? in the vicinity of /3c (right panel) point to a first order 
transition. 



In an earlier work we derived a 3 dimensional effective theory for the dynamics of the Polyakov loop for finite temperature 
G2 gluodynamics and analyzed the resulting Landau-type theory with the help of elaborate Monte Carlo simulations [14]. 
Aheady the leading order effective Polyakov loop model exhibits a rich phase structure with symmetric, ferromagnetic, and 
anti-ferromagnetic phases. 

In the present paper we investigate the phase structure of microscopic G2 YMH lattice theory with a Higgs field in the 7 
dimensional representation. The corresponding lattice action for the G2 valued link variables and a normalized Higgs field with 
7 real components reads 

5YMH[Z^,$]=/3^(l-ytrRcWn') = (4) 

and depends on the inverse gauge coupling /3 and the hopping parameter k. For /3 ^ 00 the gauge bosons decouple and the 
theory reduces to an 0(7) invariant nonlinear sigma model which is expected the have a second order (mean field) symmetry 
breaking transition down to 0(6). The mean field prediction for the critical coupling is Kc,mf — 7/8 and this value bounds Kc 
from below [15]. In the limit k = we recover G2 gluodynamics with a first order deconfinement phase transition, in agreement 
with the findings in [16]. In the other extreme case k ^ 00 we end up with SU{3) gluodynamics with a weak first order 
deconfinement transition. The known transitions in the limiting cases k — s> 0, k — > 00 or /3 — > cxd are depicted in Fig. 2. If k is 
lowered from 00 then in addition to the 8 gluons of SU{3), the 6 additional gauge bosons of G2 with decreasing mass begin to 
participate in the dynamics. Similarly as dynamical quarks and anti-quarks, they transform in the representations {3} and {3} 
of SU{3) and thus explicitly break the Z3 center symmetry. As in QCD they are expected to weaken the deconfinement phase 
transition. Thus it has been conjectured in [6] that there may exist a critical endpoint where the transition disappears. 
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In the following section we shall briefly recall those facts about G2 representations which are relevant for the present work. 
In Sec. Ill some algorithmic aspects are reviewed. A more detailed presentation can be found in our earlier paper [13]. Sec. IV 
contains our Monte-Carlo results for the phase diagram in the (/3, n) plane. We find that the two first order lines emanating from 
the deconfinement transitons in G2 and SU{i) gluodynamics at k = and n = 00 end in the vicinity of k) = (9.4, 1.6) 
on a 6 X 16'^ lattice. Sec. VI contains the results of our high statistics simulations for histograms and susceptibilities in the 
small region in parameter space where the two first order lines are either connected by a second order line or leave open a gap 
which smoothly connects the confinend and deconfined phases. Our data are consistent with the conjectured critical endpoints 
attached to the two first order lines. For large (3 a second order transition line which separates the 0(7) and 0(6) sigma models 
comes close to the first order deconfinement transition lines . The phases and transition lines are localized and analysed with high 
statistics simulations of the Polyakov loop distribution and susceptibility, plaquette and Higgs action susceptibilities, and finally 
with derivatives of the mean action with respect to the hopping parameter Besides the transition lines indicated in Fig. 2 there 
exists another line of monopole driven bulk transitions. This Une emanates from the bulk crossover in pure G2 -gluodynamics at 
/3 = 9.45 [16]. 



11. THE GROUP G2 



The exceptional Lie group G2 is the smallest Lie group in the Cartan classification which is simply connected and has a trivial 
center The two fundamental representations are the 7 dimensional defining representation {7} and the 14 dimensional adjoint 
representation {14}. One may view the elements of the representation {7} as matrices in the defining representation of SO{7), 
subject to seven independent cubic constraints, see [8]. For example, the defining representation {7} of S0{7) turns into 
an irreducible representation of G2, whereas the adjoint representation {21} of 5*0(7) branches into the two fundamental 
representation {14} and {7} of G2. The gauge group of strong interaction is a subgroup of G2 and the corresponding coset 
space is a sphere [17], 

G2/SU{3) - S'^. (5) 

This means that every element U of G2 can be written as 

U = S-V with SeG2/SU{3) and V e SU{3), (6) 

and we shall use this decomposition to speed up our numerical simulations. 

Quarks in G2 transform under the 7 dimensional fundamental representation, gluons under the 14 dimensional fundamental 
(and adjoint) representation. To better understand G2 gluodynamics we recall the decomposition of tensor products 

{7} (g) {7} = {1} e {7} ® {14} ® {27}, 
{7} {7} (g) {7} = {1} © 4 ■ {7} © 2 ■ {14} © 3 • {27} © 2 • {64} © {77'}, 

{14} (g) {14} = {1} © {14} © {27} © {77} © {77'}, ^'^^ 
{14} (g) {14} (g) {14} = {1} © {7} © 5 • {14} © 3 • {27} © ■ • ■ . 
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These decompositions show similarlies to QCD: two quarks, three quarks, two gluons and three gluons can build colour singlets 
- mesons, baryons and glueballs. In G2 gauge theory three gluons can screen the colour charge of a single quark, 

{7}«){14}(»{14}®{14} = {1}®--- , (8) 

and this explains why the string between two external charges in the {7} representation will break for large charge separations. 
The two remnants are colour blind glue lumps. The same happens for two external charges in the adjoint representation. In a 
previous work we did observe string breaking at the expected separation between the two charges [13]. 

The G2 gauge symmetry can be broken to 5?7(3) with the help of a Higgs field in the 7 dimensional representation. For 
K ^ 00 the factor S in the decomposition (6) is frozen and we end up with an SU{'i) gauge theory with rescaled gauge coupling 
for the factor U. With respect to the unbroken subgroup SU{i) the fundamental representations {7} and {14} branch into the 
following irreducible SU{3) representations: 

{7}^{3}®{3}®{1}, 

{i4}^{8}e{3}e{3}. 

The Higgs field branches into a scalar quark, scalar anti-quark and singlet with respect to SU{3). Similarly, the G2 gluons 
branch into massless SU{3) gluons and additional gauge bosons with respect to SU (3). The latter eat up the non-singlet scalar 
fields such that the spectrum in the broken phase consists of 8 massless gluons, 6 massive gauge bosons and one massive Higgs 
particle. 



in. ALGORITHMIC CONSIDERATIONS 
A. Equations of motion for local hybrid Monte-Cario 



In this work we employ a local version of the hybrid Monte-Carlo (HMC) algorithm where single site and link variables are 
evolved in a HMC style [18]. The algorithm assumes a local interaction and hence applies to all purely bosonic theories. The 
implementation for the G2 Higgs model is a mild generalization of the algorithm used in our previous work on G2 gluodynamics 
[13]. We use a local hybrid Monte-Carlo (LHMC) algorithm for several good reasons: First there is no low Metropolis acceptance 
rate even for large hopping parameters. More precisely, in a heat bath algorithm combined with an over-relaxation we would 
need two Metropolis steps in each update for k > which for large k may lead to low acceptance rates. With the LHMC- 
algorithm we can avoid this problem and deal with arbitrary values of k. Autocorrelation times can be controlled (in certain 
ranges) by the integration time in the molecular dynamics part of the HMC algorithm. Second, the formulation is given entirely 
in terms of Lie group and Lie algebra elements and there is no need to back-project onto the group. For G2 it is possible to use 
a real representation and in addition an analytical expression for the involved exponential maps from the algebra to the group. 
These maps allow for a fast implementation of the LHMC algorithm. 

This algorithm has been essential for obtaining the accurate results in the present work. Since we developed and used the first 
implementation for G2 it may be useful to sketch how it works for this exceptional group. More details can be found in [13]. 
For G2 YMH lattice theory the (L)HMC algorithm is based on a fictitious dynamics for the link-variables on the G2 manifold 
and the normalized Higgs field on the 6-sphere. The "free evolution" on a semisimple group is the Riemannian geodesic motion 
with respect to the Cartan-Killing metric 

ds% = Ktr {dUU-^ (g) dUU-^) . (10) 

In a (L)HMC dynamics the interaction term is given by the YMH action (4) of the underlying lattice gauge theory and hence it 
is natural to derive the HMC dynamics from a Lagrangian of the form 

Lhmc = (^x.M^^.-i)' + K{^,^) - Symh[U,^], (11) 

where 'dot' denotes the derivative with respect to the fictitious time parameter r and K(^, $) is a kinetic term for the Higgs 
field. To update the normalized Higgs field we set 

= O^$o with e SO{7) (12) 

and constant $o- The change of variables Ox converts the induced measure on C into the Haar measure of SO {7). 

Without interaction the rotation matrices Ox will evolve freely on the group manifold 50(7) such that in terms of the {U, O) 
variables we choose as Lagrangian for the HMC dynamics 

L = -\Y.^^ {i^x,^.U-iy -l^tv (OxO-'Y - 5ymh[W, O] . (13) 

X,fJ. X 
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The Lie algebra valued fictitious momenta conjugated to the link variable Ux,ij. and site variable Ox are given by 

The Legendre transform yields the following pseudo-Hamiltonian 

H = -1 ^trq}2_^ _ i ^trn^ + 5ymh[Z^, O]. (15) 

Note that for real U^,^ and the momenta are antisymmetric such that both kinetic terms are positive. The equations of motion 

for the momenta are obtained by varying the Hamiltonian. The variation of 5ymh[ O] with respect to a fixed link variable Ux.^ 
yields the staple variable Rx,^, the sum of triple products of elementary link variables closing to a plaquette with the chosen link 
variable. Setting 

S^x,^l = ^x,fj.dT, 5Ux,n = Ux,t_idT = -^x,piUx,fi,dT (16) 

with similar expressions for the momentum and field variables Q.x and Ox in the Higgs sector yields for the variation of the 
HMC Hamiltonian 

5H = - ^trq3,,4*P,,^ - Fx,^] - ^trO^Q, - G,} (17) 

X^IJ, X 

with the following "forces" in the gauge and Higgs sector 

Fx,l. = ^ {Ux,^.Rx,^. - Rl^Uly) + KiUxM<Pl+^, Gx = ^<Px(j2^,^l(xy <^j;)^, (18) 

where the last sum extends over all nearest neighbors y of x and Uxy denotes the parallel transporter from y to x. The variational 
principle implies that the projection of the terms between curly brackets onto the Lie algebras 52 and so (7) vanish, 

^x,n = Fij.,x\q^ j ~ ^^\so{7)' ^^^^ 

The equations (14) and (19) determine the fictitious dynamics of the lattice fields in the (L)HMC algorithm. Choosing a trace- 
orthonormal basis {Ta} of Q2 the LHMC equations in the gauge sector read 

Ux,^. = -'^x,^Mx,^. and q?x,^ = ^ tr (F^,^T„) T„ (20) 

a 

with force Fx^^ defined in (18). In the Higgs sector they take the form 

dx = -i2xOx and Qx = Y.^r(Gxfb^n (21) 

6 

with trace-orthonormal basis {T^} of so(7) and force Gx defined in (18). 



B. Numerical solutions of YMH-dynamics 

We employ a time reversible leap frog integrator which uses the integration scheme 

"^xA^ + s^) = "^^A^ + k^^) + h^^VxA^ + ^r), 

and similarly for the variables {Ox,Q.x) in the Higgs sector. The 'time' derivative of *P(t + St) in the last step is given in 
terms of the already known group valued field at t + St via the equations of motion. Clearly, to calculate U and O at time 
T + St a fast implementation of exponential maps is required. In the Higgs sector the map so(7) — > 50(7) is computed via the 
Cayley-Hamilton theorem. For small values of the hopping parameter k the step size and integration length for the integration 
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may be chosen as in the gauge field integrator. For an efficient and fast computation of the exponential map 02 G2 we exploit 
the real embedding V of the representation 3 © 3 of SU{3) into G2, 

U = S-V{W) with SgG2/SU{3), W e SU{3). (23) 

For a given time step St the factorization will be expressed in terms of the Lie algebra elements with the help of the exponential 

maps, 

exp {St u} = exp {6ts} ■ exp {6t 0} with generators u € 32, t) e V* (su(3)). (24) 

The exponential maps for the two factors can be calculated efficiently, see [13]. But in the numerical integration we need the 
exponential map for elements u € 02- These elements are related to the generators s and D used in the factorization by the 
Baker-Campbell-Hausdorff formula, 

6tu = St{s + \:>) + ^St^[s,v] + --- . (25) 

For a second order integrator the approximation (25) may be used in the exponentiations needed to calculate V and S. This 
approximation leads to a violation of energy conservation which is of the same order as the violation one finds with a second 
order integrator. To sum up, a LHMC sweep consists of the following steps: 

1. Gaussian draw for the momentum variables on a given site and link, 

2. Integration of the equations of motion for the given site and link, 

3. Metropolis accept/reject step, 

4. Repeat these steps for all sites and links of the lattice. 

This local version of the HMC does not suffer from an extensive 6H cx V problem such that already a second order symplectic 
(leap frog) integrator allows for sufficiently large time steps St. For a large range of couplings k) in our simulations an 
integration length of T = 0.75 with a step size of St = 0.25 is optimal for minimal autocorrelation times and a small number of 
thermalisation sweeps. Acceptance rates of more than 99% are reached. To compare the performances of our LHMC algorithm 
with the usually used heat-bath algorithm we estimated the computation time of the different parts in the LHMC-algorithm in 
units given by the average computation time for one staple in ASu- On an Intel Corei? CPU the latter is approximately 4 fis for 
a 12^ X 6 lattice. 

In Table I we listed the times needed to change the gauge or Higgs action during a single update of one link or one Higgs field 
variable, the time for both integrators without exponential map and separately the computation time for a single exponential map. 
Most time is spent with calculating the exponential maps for 50(7). Note that during the calculation of one exponential map for 

Part ASu ASp integr. U integr. O exp(G2) exp(5'07) 

pure gauge 1.00 - 1.34 - 0.42 

gauge Higgs 1.03 0.43 1.74 1.00 0.40 4.97 

TABLE I: Computation times normalized to ASu (staple). 

SO{7) the CPU calculates about 10 exponential maps for G2. Table II compares the total time-contributions to one configuration 
with those of the heat-bath algorithm with overrelaxation. We see that for pure gauge theories the standard heat-bath algorithm 
with overrelaxation is only two times faster as the LHMC algorithm. 

Part ASu ASp integr. U integr. O exp(G2) exp(g07) total time/V ■ d ■ Config. heat-bath 

pure gauge 1.00 - 1.34 - 1.26 - 3.60 w2 

gauge Higgs 1.03 0.11 1.74 0.25 1.20 3.72 8.05 

TABLE II: Total time contribution to one LHMC configuration compared to heat-bath algorithm. 
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IV. THE PHASE DIAGRAM OF THE G2 HIGGS MODEL: OVERVIEW 

With the help of the local HMC algorithm sketched previously we calculated several relevant observables to probe the phases 
and phase transition lines in the n) plane. First we present the phase diagram obtained on small lattices. For vanishing k we 
are dealing with G2 gluodynamics which shows a first order finite temperature deconfinement phase transition. The transition 
is discontinuous since there is a large mismatch of degrees of freedom in the confined and unconfined phases. At the other 
extreme value k = cxd six of the fourteen gauge bosons decouple from the dynamics and we are left with SU (3) gluodynamics, 
which shows a first order deconfinement phase transition as well. The question arises whether the first order transitions in G2 
and SU (3) gluodynamics are connected by a unbroken line of first order transitions or whether there are two critical endpoints. 
In the latter case the confined and unconfined phases could be connected continuously. On the other hand, for arbitrary n but 
/3 — > 00 the gauge degrees of freedom decouple from the dynamics and one is left with a nonlinear 0(7)-sigma model. We 
expect that the 0(7)-symmetry is spontaneously broken to 0(6) for sufficiently large values of the hopping parameter and that 
this transition is of second order 

In order to localize the confinement-deconfinement transition line(s) we first measured the Polyakov loop expectation value 
as (approximate) order parameter for confinement on a small 12"^ x 2-lattice in a large region of parameter space (/3 = 5 ... 10, 
K = . . . 10'*). For K 1 the Polyakov loop takes its values in the reducible representation {3} © {3} © {1} of 5?7(3) and 

{P)^l + {P + P)su(z)- (26) 

Thus, for large k we should find {P) « 1 in the confining phase and {P) « 7 or {P) w —2 in the unconfined phase where P 
is near one of the three center-element of SU (3). We eliminate the ambiguity of assigning a value to the Polyakov loop in the 
unconfined phase by mapping values with (P) < 1 to 3 — 2 {P). 

The result for {P) is depicted in Fig. 3. We see that in the confining phase the expectation value varies from to 1 when 




the hopping parameter increases. For large values of /3 in the unconfined phase the Polyakov loop is near the identity or (for 
large k) near one of the three center-elements of SU (3). On the small lattice the Polyakov loop jumps along a continuous curve 
connecting the confinement-deconfinement transitions of pure G2 and pure SU (3) gluodynamics. This suggests that there exists 
a connected first order transition curve all the way from k to k = cx). To see whether this is indeed the case we performed 
high-precision simulations on larger lattices. A careful analysis of histograms and susceptibilities for Polyakov loops and the 
Higgs action shows that the first order lines beginning at k = and at k — cx) do not meet. This happens in a rather small region 
in parameter space such that the two first order lines almost meet. They may be connected by a line of continuous transitions or 
in-between there may exists a window connecting the confined and unconfined phases smoothly. 
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For /? — cxD we are left with a nonlinear 0(7) sigma model with action 

S„ = -nJ2'^^+i^^^^ (27) 

and this model shows a second order transition at a critical coupling Kc from a 0(7) symmetric to a 0(6) symmetric phase. To 
see how this transition continues to finite values of f3 we measured the expectation values (Op) and (Oh) of the (averaged) 
plaquette variable and Higgs action 

Op = RctrWn and Oh = • (28) 

□ XfJ, 

and the corresponding susceptibihties 

xiO) = V{{0')-{Of) . (29) 

The finite size scaling theory predicts that near the transition point the maximum of the susceptibilities scales with the volume 
to the power of the corresponding critical exponent 7 

xiO) - ar/'' + b , (30) 

where i' is the critical exponent related to the divergence of the correlation length. For a first order phase transition we expect 
the susceptibility peak to scale linearly with the spatial volume (since Nf is fixed). More precisely, for a first order transition one 
expects 7 = 1 and = 1/3 while for a second order transition 7 7^ 1 [19]. 

The expectation values and logarithms of susceptibilities on a small 6^ x 2-lattice are depicted in Fig. 4. The expectation 




FIG. 4: Average plaquette, Higgs action and susceptibilities near the critical point on 6'^ x 2 lattice. 



value of a plaquette variable jumps at the deconfinement transition line and the corresponding susceptibility is peaked. This 
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is in full agreement with the jump of the Polyakov loop across this transition line. The expectation value of the Higgs action 
and the corresponding susceptibility both spot the deconfinement transition well. But they also discriminate between the 0(7) 
unbroken and broken phases. The data on the small lattice point to a second order Higgs transition line in the YMH-model for all 
P > /^dcconf ('*)■ This could imply that the second order line ends at the first order deconfinement transition line. To determine 
the order of the Higgs transition line we consider the finite size scaling of 

x{OH)^-g^{OH) and q^{Oh) (31) 

for lattices up to 20'^ x 6. The results presented below show that the Higgs transitions are second order transitions. Unfortunately 
we cannot exclude the possibility that the second order line turns into a crossover near the deconfinement transition line. 

Our results on the complete phase diagram in the {13, K)-plane as calculated on a larger 16'^ x 6-lattice are summarized in 
Fig. 5. We calculated histograms and susceptibilities near the marked points on the transition lines in this figure. If the triple 
point exists then an extrapolation to the point where the confined phase meets both unconfined phases leads to the couplings 
/3(jip = 9.62(1) and Ki^p — 1.455(5). Near this point the deconfinement transition is very weak, continuous or absent and thus 

oo I m» 1 1 1 n 2 I i«H [-» 1 1 1 



so {7) broken 
SU (3) deconfinement 



50(7) unbroken 
G2 deconfinement 



10 15 
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FIG. 5: Phase transition lines on a 16 x 6 lattice. The solid line corresponds to the first order deconfinement transitions, the dashed line to 
the second order Higgs transitions and the dotted line to the left of the first order line to the bulk transitions. The plot on the right shows the 
details inside the small box in the plot on the left where the transition lines almot meet. 



we performed high-statistics simulations on larger lattices to investigate this region in parameter space more carefully. Some 
of our results are presented in the following sections. Up to a rather small region surrounding (/3trip7 '^tnp) we can show that the 
deconfinement transition is first order and the Higgs transition is second order But we shall see that in a small region around 
this point the deconfinement transition is either second order or absent. 



The bulk transition 



The existence of a bulk transition in lattice gauge theoiies at zero temperature can influence its finite temperature behaviour 
Such transitions are almost independent of the size of the lattice and are driven by lattice artifacts [20]. Bulk transitions between 
the unphysical strong-coupling and the physical weak-coupling regimes in lattice gauge theories is the rule rather than the 
exception. The strong coupling bulk phase contains vortices and monopoles which disorder Wilson loops down to the ultraviolet 
length scale given by a^a ^ 0{1) [21, 22]. In the weak coupling phase the short distance physics is determined by aymptotic 
freedom and a^a -C 1. Both SU{2) and SU (3) lattice theories exhibit a rapid crossover between the two phases which beomes 
more pronounced for SU{4) [21]. For SU{N) with N > 5 the bulk transition is first order [21]. SU{3) lattice gauge theory 
with mixed fundamental (/) and adjoint (a) actions shows a first order bulk transiton for large /3a and small For decreasing 
(3a the transition line terminates at a critical point and turns into a crossover touching the line (3a = 0. On lattices with Nt = 2 
the deconfinement transition line joins the bulk transition line smoothly from below and for Nt > 4: from above [23, 24]. More 
relevant for us is the finding in [16] that the bulk transition in pure G2 gauge theory at (3 ~ 9.45 is a crossover [16]. 

We have scanned the values for the plaquette variables and Polyakov loops from the strong to the weak coupling regime to 
find a bulk transition that might interfere with the finite temperature deconfinement transition. For various values between k = 
and K = cx) on a 12'^ x 6 and 16'^ x 6 lattice we determined the position and nature of the bulk transitions. In full agreement 
with [16] we see a crossover at (/3, k) « (9.44, 0) which is visible as a broad peak in the plaquette susceptibility depicted in the 
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FIG. 6: Plaquette and susceptibility for small values of k near the bulk transition on a 12'^ x 6 lattice. 



right panel of Fig. 6. The Polyakov loop does not detect this crossover. Note that for small k the position of the bulk transition 
does not depend on the hopping parameter which means that the bulk transition line hits the line k ~ vertically. Despite of 
the broad peak in the susceptibility of the plaquette density are the bulk and deconfinement transition cleary separated and this 
agrees with the results in [25]. In the region 1.3 < k < 1.6 the critical coupling f3c decreases with increasing Kc but the nature 
of the transition does not change much as can bee seen in Fig. 7. The plaquette density seems to be a continous function of (3 
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FIG. 7: Plaquette and susceptibility for intermediate values of k near the bulk transition on a 12^ X 6 lattice. 



and K and we conclude that the transition is still a crossover 

Between k = 1.6 and k — 1.65 the peak in the bulk transition becomes pronounced. In this region the distance between 
the bulk and deconfinement transitions becomes very small. Nevertheless we expect that the much localized bulk transition still 
does not interfere with the weak deconfinement transition. For values of k between 1.65 and approximately 2.5 the position of 
the bulk transition gets more sensitive to the hopping parameter and the distance to the deconfinement transition line increases 
again. The nature of the transition changes at the same time - a large gap in the action density separates the strong coupling 
from the weak coupling region. This is depicted in Fig. 8. The many data points taken at k = 1.8 show that the size of the gap 
does not depend on the volume and this points to a first order transition. The plots for the plaquettes and plaquette susceptibilites 
look very much like the plots in Fig. 6. For k > 2.5 the situation changes again. The gap in the plaquette density closes and the 
position of the bulk transition tends to that of the bulk transition in SU (3) gluodynamics which again is a crossover 

There is ample evidence that bulk transitions are driven by monopoles on the lattice [20]. Thus we calculated the density of 
monopoles [24] as a function of (3 for k = and k = 1.8. The density M together with the plaquette variable are plotted in 
Fig. 9. For k = they vary smoothly with (3, as expected for a cross-over, but for k = 1.8 they jump at the same /3 w 9.25. 
The height of the jump does not depend on the lattice size, see Fig. 9, right panel. Thus we find strong evidence that the bulk 
transition is intimately related to the condensation of monopoles in the strong coupling G2 Higgs model. 

Finally we would like to comment on the behaviour near k — 1.6. Here the G2 Higgs model behaves similar to SU{3) 
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gluodynamics with mixed fundamental and adjoint actions. The latter shows a first order bulk transition which turns into a 
crossover for small fia- It seems that for n > 1.6 the massive G'2-gluons are heavy enough such that the approximate center 
symmetry of the unbroken SU{3) is at work. This could explain why we find a first order transition for k> 1.6. 



V. THE TRANSITION LINES AWAY FROM THE TRIPLE POINT 



In this section we come back to the confinement-deconfinement transition. Sufficiently far away from the suspected triple 
point at /Strip = 9.62(1) and Ktdp = 1.455(5) the signals for first- and second order phase transitions are unambiguous and are 
presented in this section. The measurements taken near the would-be triple point are less conclusive and will be presented and 
analysed in the following section. 



The conlinement-deconflnement transition line 



Already the histograms for the Polyakov loop show that the deconfinement transition is first order for values of the hopping 
parameter k in the intervals [0, 1.4] and [1.7, 00]. Two typical distributions for k — 1.0 and k = 1.3 corresponding to the points 
1 and 2 in the phase diagram in Fig. 5 are depicted in Fig. 10 (left panel). These and other histograms with k < 1.4 show a clear 
double peak structure near the transition line and are almost identical to the histogram for k = 0. Similar results are obtained 
for larger hopping parameters k > 1.7. 

In Fig. 10 (Right panel) we plotted histograms of the Polyakov loops for /3 = 9 and hopping parameters in the vicinity of 
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FIG. 10: Distributions of the Polyakov loop on a 16^ x 6-lattice. Left panel: {/3, k) = (9.76, 1) and (9.725, 1.3); Right panel: /3 = 9 and 
various values of the hopping parameter. 

K « 2.6, corresponding to point 3 in Fig. 5. The histograms with k < 2.6525 show peaks at almost the same positions. The 
systems with these small values of k are in the confined phase. For larger K-values the peak moves towards the 'would-be' center 
elements of the subgroup SU{3) and a second peak appears. Again the double-peak structure of the distribution points to a first 
order transition. We varied the spatial sizes of the lattices and observed no finite size effects in the distributions for Ng > 16. 



The Higgs transition line 



For P ^ oo the gauge degrees of freedom are frozen and we are left with a nonlinear 0{7) sigma-model which shows a 
second order transition from a 0(7)-symmetric massive phase to a 0(6)-symmetric massless phase. With the help of a cluster 
algorithm [26] we updated the constrained scalar fields and calculated the susceptibility of 



which is proportional to the sigma-model action in (27), 

x(a) -- 



XjfJ, 



kV 



(32) 



(33) 



The results of our simulations on lattices with varying spatial sizes are depicted in Fig. 11, left panel. The susceptibility of 
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FIG. 11: The first and second derivative of the average sigma-model action for different spatial lattice sizes. 
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the action becomes steeper as the spatial volume increases while the peak of the (normalized) second derivative also increases. 
This means that the system undergoes a second order transition at — 1.075(5) (corresponding to point 4 in Fig. 5) from 
a massive 0(7)-symmetric phase with vanishing vacuum expectation value to a massless 0(6)-symmetric phase with non- 
vanishing expectation value. Actually the mean field theory for 0{n) models in d dimensions predicts a second order transition 
at the critical coupling Kc.mf ~ n/2d. For our model in 4 dimensions the mean-field prediction is Kc.mf = 7/8 « 0.875 and is 
not far from our numerical value. 

For smaller values of f3 the gauge degrees of freedom participate in the dynamics and 3^ (S) is now proportional to the 
susceptibility of Oh in (28). The plots in Figs. 12 and 13 show a similar behavior of the first and second derivatives of the 
average Higgs action for /3 — 30 and 12, corresponding to the points 5 and 6 in the phase diagram in Fig. 5. 
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FIG. 12: First and second derivative of the average action with respect to the hopping parameter for different spatial lattice sizes at /3 = 30 
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FIG. 13: First and second derivative of the average action with respect to the hopping parameter for different spatial lattice sizes at /3 = 12 



Even for the smaller value /3 = 12 we see that the susceptibility becomes steeper with increasing lattice size while the second 
derivative of the average action increases. This already demonstrates that the second order transition at the aymptotic region 
/3 — cx) extends to smaller values of /3. 



VI. THE TRANSITION LINES NEAR THE TRIPLE POINT 



When the first order transition become weaker it becomes increasingly difficult to distinguish it from a second order transition 
or a cross-over For example, the four histograms in Fig. 14 show distributions of the Polyakov loop at point 7 in the phase 
diagram depicted in Fig. 5, corresponding to k = 1.5 and /3 varying between 9.5525 and 9.5550. All histograms are computed 
from 400 000 configurations on a medium size 16'^ x 6 lattice. The histogram on top left shows a pronounced peak at P « 0. 1, 
corresponding to the value in the confined phase. With increasing f3 a second peak builds up at P « 0.25 corresponding to a 
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value in the unconfined phase. We have calculated more histograms and conclude that the well-separated peaks in the distribution 
are of equal heights for (3^ ~ 9.5535. At this point the Polyakov loop jumps from the smaller to the larger value. For even larger 
values of (3 the second peak at larger P takes over and the system is in the unconfined phase. Although the histograms point to a 




FIG. 14: Distributions of the Polyakov loop at «: = 1.5 where the transition is weakly first order on a 16^ x 6 lattice with 400 000 configurations 
for each histogram. Top left /3 = 9.5525, top right /3 = 9.5535, bottom left /3 = 9.5540 and bottom right /3 = 9.5550 (JS^ ^ 9.5535). 

weakly first order transition we can not rule out the possibility that the transition at k = 1.5 and /3 w 9.5535 is of second order 
Later we shall see that it is a first order transition. If we slightly decrease the value of k, then the signal for a first order transition 
is more pronounced. This is illustrated in the Polyakov loop histograms depicted in Fig. 15. If we again increase the value from 
K = 1.5 to K = 1.55 the peak of the Polyakov loop does not jump at the transition point at /3 w 9.4885. Instead it increases 
smoothly from P « 0.12 in the confinement phase to P w 0.24 in the deconfinement phase, see Fig. 16. We conjecture that in 
this region of parameter space the first order transition turns into a continuous transition or a cross-over which is later confirmed 
by an even more careful analysis. 

We studied the size-dependence of the average Polyakov loop, plaquette variable and Higgs-action per lattice site together 
with their susceptibilities. The following results are obtained on lattices with Nt = 6 and spatial extends Ng G {12, 16, 20, 24} 
and for /? = 9.5535. This corresponds to points in the neighborhood of point 7 in the phase diagram in Fig. 5. 

Fig. 17 shows the K-dependence of the Polyakov loop and its susceptibility for the four different lattices. The measurements 
have been taken at 20 different values of the hopping parameter in the vicinity of k = 1.5. This way we cross the phase transition 
line vertically in the K-direction at the transition point 7 in the phase diagram in Fig. 5. The K-dependence has been calculated 
with the reweighting method. Later we shall see that the peak of the susceptibility at ~ 1.501 scales linearly with the volume. 
This linear dependence is characteristic for a first order transition. 

The plots in Fig. 18 show the K-dependence of the average plaquette variable and the corresponding susceptibility for the four 
lattices. Again we observe that the susceptibility peak at Kc « 1.501 increases linearly with the volume of the lattice. Also note 
that on the small 12^^ x 6 lattice the peak in the susceptibility can hardly be seen. 

The two plots in Fig. 19 show the K-dependence of the average Higgs action per lattice point and corresponding susceptibility. 
Similarly as for the Polyakov loop and the plaquette we observe a peak of the susceptibility at the same value Kc ~ 1.501. 

To check for finite size scaling we investigated the susceptibilities corresponding to the Polyakov loop, plaquette variable and 
Higgs-action per site as a function of the volume. The results are plotted in Fig. 20, left panel. For an easier comparison we 
normalized the data points by the peak value for the largest lattice with lattice size Ng — 24. The linear dependence of the peak 
susceptibilities on the volume is clearly visible for the larger three lattices and this linear dependence is predicted by a first order 
transition [19]. In recent studies of the lattice SU (2) Higgs model in [25] it turned out that for Ng = Nt < 18 the maxima of the 
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FIG. 15: Distribution of the Polyalcov loop at (/3, k) = (9.6190, 1.455) - (9.6220, 1.455) near the supposed triple point; 400 000 configura- 
tions on 16'' X 6 lattice 

susceptibilities are well described by a function of the form aL'^ + b, so that they seem to scale linearly with volume, as expected 
for a first order transition at zero temperature. Simulations on larger lattices revealed however, that the suceptibility peaks all 
saturate at larger values of L and no singularities seems to develop in the thermodynamic limit. For the lattice G'2-Higgs model 
considered in the present work we see no flattening of the peaks for larger lattices with Ns up to 24 and we interpret this as a 
signal for a true first order transition. 

Table III shows the extrapolation of the critical hopping parameter to infinite volumes. To that end we calculated for each 
lattice size the value Kc{V) at which the Polyakov loop-, plaquette- and Higgs action susceptibilities take their maxima. Note 
that on the larger lattices with Ns = 20 and 24 the three critical hopping parameters are the same within statistical errors. The 
infinite volume extrapolation yields the critical value Kc = 1.5008. 

Volume 12^^ 16^ 20^ 24^ 
xiP) 1.5012 1.5016 1.5008 1.5008 
xiOn) 1.4992 1.5012 1.5008 1.5008 
xiOp) 1.4980 1.5008 1.5008 1.5008 

TABLE III: Critical coupling tic obtained from the maximum of the susceptibility peaks of Polyakov loop, plaquette and Higgs action for 
different spatial volumes al P — 9.5535 



A. The first order lines do not meet 

The previous results on the 16^ x 6 lattice leave a small region in parameter space near (/3, k) (9.4, 1.6), where the transition 
may be continuous or where we can cross smoothly between the confined and unconfined phases. Since a jump of the Polyakov 
loop expectation values in the infinite volume limit points to a first order transition we investigated the quantity 
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FIG. 16: Distributions of the Polyakov loop at «: = 1.55 where the transition is probably not first order on a 16 x 6 lattice with 800 000 



configurations for each histogram. Top left /? — 9.4875, top right /? 
(Pc ~ 9.4885). 



= 9.4885, bottom left /3 = 9.4895 and bottom right /3 = 9.4905 
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FIG. 17: Finite size scaling of Polyakov loop and Polyakov loop susceptibility at /? = 9.5535 



more carefully. In the small parameter region we localized the critical curve {(3c, Hc) with the histogram method. At the critical 
point is the height of the confinement peak equal to the height of the deconfinement peak. For fixed we crossed the transition 
line by increasing the inverse gauge coupling. Then we measured the maximal jump as a function of the step size A/3 for one 
step size below and one above f3c. For a first order transition the jump should not depend much on A/3 whereas for a continuous 
transition or a cross-over AP should decrease with decreasing A/3. The results on a 16'^ x 6 lattic are depicted in Fig. 20 (right 
panel). We see that for 9.35 < /3c < 9.52 corresponding to 1.52 < < 1.72 the jump approaches zero with shiinking step size 
and this clearly points to second order confinement-deconfinement transitions or cross-overs in these small parameter regions. 
Simulations on a larger 20^ x 6 lattice confirm these results. Fig. 21 shows histograms of the Polyakov loop for K-values between 
1.5 and 1.7. At k = 1.5 we still observe a weakly first order transition which turns into a continuous transition or crossover 
for 1.5 < K < 1.7. Within the given resolution in parameter space the window is the same as on the 16^ x 6 lattice. Since the 
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FIG. 18: Finite size scaling of the plaquette variable and its susceptibility for /3 = 9.5535. 
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FIG. 19: Finite size scaling of Higgs action and its susceptibility for /3 = 9.5535. 



critical couplings for spatial volumes beyond 20'^ do not change we conclude that the gap will not close in the infinite volume 
limit. This shows that the two first order lines emanating from k = and k = oo do not meet. 

Here the question arises whether such a gap in the first order line between the confined and unconfined phases is expected. The 
celebrated Fradkin-Shenker-Osterwalder-Seiler theorem [27, 28], originally proven for the SU{N) Higgs-model with scalars in 
the fundamental representation, says that there is no complete separation between the Higgs- and the confinement regions. Any 
point deep in the confinement regime and any point deep in the Higgs regime are related by a path such that Green's functions 
of local, gauge invariant operators vary analytically along the path. Thus there is no abrupt change from a colorless to a color- 
charged spectrum. This is consistent with the fact that there are only color singlet asymptotic states in both 'phases'. 

The proof of the theorem relies crucially on using a completely-fixed unitary gauge. A complete gauge fixing is not possible 
with scalars in the adjoint representation of SU (N) since these scalars are center blind. Thus the theorem does not hold for 
adjoint scalars and indeed, with adjoint scalars there exits a phase boundary separating the Higgs and confined phases. It is not 
completely obvious what these results tell us about the phase diagram of the G2 Higgs model. The center of G2 is trivial and 
the 14-dimensional adjoint representation is just one of the two fundamental representations. Since there is no need to break the 
center one may conclude that the confinement-like regime and the Higgs-like regimes are analytically connected. In addition, 
for large values of the hopping parameter the center of the corresponding SU{3) gauge theory is explicitly broken by the scalar 
fields, simililarly as for the SU{3) Higgs model with scalars in the fundamental representation. These arguments suggest that 
there exist a smooth cross-over between the confining and Higgs phases. But one important assumption of the Fradkin-Shenker 
theorem is not fulfilled for the G2 Higgs model. The theorem assumes that there exists no transition for large k. Then at large k 
one can move from large to small f3 and then at small f3 further on to small values of k without hitting a phase transition. Clearly 
this is not possible for the G2 Higgs model such that not all assumption of the theorem hold true. 
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FIG. 20: Left panel: finite size scaling of the three susceptibilities at the transition point with /3 — 9.5535. The lines are fits to the peak values, 
Xmax(V^) = aV + b. Right panel: Difference of Polyakov loop in confined and unconfined phase at the phase transition point for various 
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FIG. 21: Distribution of Polyakov loop near the phase transition point for k = 1.5 (top left), k = 1.55 (top right), k = 1.65 (bottom left) and 
K = 1.7 (bottom right) on a 20^ x 6 lattice 



19 



VII. CONCLUSIONS 



With a new and fast LHMC-implementation for the exceptional G2 Higgs model we calculated the full phase diagram in the 
coupling constant plane spanned by the hopping parameter n and inverse gauge coupling /3. First we confirmed the proposed 
and earlier seen [6, 16] first order transition for pure G'2-gluodynamics which corresponds to the line k = in the phase 
diagram of the Higgs model. A first analysis on smaller lattices indicated that this first order transition is connected to the first 
order deconfinement transition in S'C/(3)-gluodynamics, corresponding to the limit k — s> 00, by a smooth curve of first order 
transitions. The same analysis spotted another curve of second-order transitions emanating from /3 — > 00 and meeting the first 
order line at a triple point. For this first analysis we calculated histograms for the Polyakov loop, Higgs-action and plaquette 
action. To identify the second order transition line we studied the finite size scaling of various susceptibilities and the second 
derivative of the action with respect to the hopping parameter. The final result of our analysis on a IQ^ x 6 lattice is depicted 
in Fig. 22. Note that the tiny region in the vicinity of the would-be triple point is very much enlarged in this figure. In this tiny 
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FIG. 22: Complete phase diagram in the (/3, «:)-plane on a 16'^ x 6 lattice. The neighbourhood of the 'would-be triple point' is very much 
enlarged and the variable scale in the diagram is responsible for the cusps in the transition lines. The solid line indicates a first order transition, 
the dashed line (blue) a second order transition and the dotted line (red) a second oder transition or a crossover. 



region in the (/?, K)-plane where the order of the transition could not be decided we studied the slope of (P) in the vicinity of 
the suspected transition. The simulations show that the two first-order curves emanating from the lines with k = and k = 00 
end before they meet. The two curves could be connected by a line of second-order transition or they could end at two (critical) 
endpoints in which case the confined and unconfined phases are smoothly connected. If indeed there exists a cross-over in G2 
Higgs model at a finite value of the hopping parameter then the gauge model behaves very similar to QCD with massive quarks. 

To finally answer the question about the behavior of G2 Higgs model theory in the vicinity of the 'would-be triple point' 
at (/?, n) K, (9.4, 1.6) further simulations with an even higher statistics and a more sophisticated analysis of the action sus- 
ceptibilities may be necessary. Since we already used an efficient (and parallelized) LHMC-algorithm and much CPU-time to 
arrive at the results presented in the work this will not be an easy task. Earlier studies of the susceptibility peaks in the simpler 
S'C/(2)-Higgs model on smaller lattices pointed to a first order transition at /? < 2.5. Recent simulations on larger lattices in [25] 
showed that the susceptibility peaks do not scale with the volume such that there is actually no first order transition for these 
small values of /?. We have seen no flattening of the peaks with the increasing volumes for Ng < 24 and conclude that the solid 
line in Fig. 22 is a first order line. But of course we cannot exclude the possibility that the correlation length is larger as expected 
and that simulations on even larger lattices are necessary to finally settle the question about the position and size of the window 
connecting the confined with the unconfined phase. This will not be easy and thus it would be very helpful to actually prove that 
the confining and Higgs phases of G2 can be connected analytically, perhaps with similar arguments as they apply to SU [N) 
Higgs models with matter in the fundamental representations [27, 28]. 
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